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ABSTRACT 

We investigate the stability of incompressible, exact, non-ideal magnetorotational (MRI) modes against par- 
asitic instabilities. Both Kelvin-Helmholtz and tearing-mode parasitic instabilities may occur in the dissipative 
regimes accessible to current numerical simulations. We suppose that a primary MRI mode saturates at an 
amplitude such that its fastest parasite has a growth rate comparable to its own. The predicted alpha parameter 
then depends critically on whether the fastest primary and parasitic modes fit within the computational domain 
and whether non-axisymmetric parasitic modes are allowed. Hence even simulations that resolve viscous and 
resistive scales may not saturate properly unless the numerical domain is large enough to allow the free evo- 
lution of both MRI and parasitic modes. To minimally satisfy these requirements in simulations with vertical 
background fields, the vertical extent of the domain should accommodate the fastest growing MRI mode while 
the radial and azimuthal extents must be twice as large. The fastest parasites have horizontal wavelengths 
roughly twice as long as the vertical wavelength of the primary. 

Subject headings: accretion, accretion disks — black hole physics — instabilities — MHD — turbulence 



1. EVITRODUCTION 

Understanding the processes that halt the ex ponential 
growth of the magnetorotat i onal instability (M RI; Velikhov 
Il95l IChandrasekh^[T96(it iBalbus & Hawlevl [T991.. 1998.) 
and set the rate of angular momentum transport in the 
turbulent regime has been an outstanding problem in ac- 
cretio n physics for almost t wo decade s . The net magnetic 
flux (iHawley etalJ 119951: ISano et alJ 2004; Pessah et al. 
2007h. the geome try of the domain (Hawlev et al. 1995; 
BodoetalJ l2008h, the re solution (Pessah et al. 2007; 
Fromang & Papaloizou' '2007*) , and the microphysical dis- 



sipati on co efficients (Flemin g et al.l [200 0; Sano & Inutsuka 
20011; lFk) mang et al. 2007; lLesur & Longarettil 120071: 
Masada & Sano 200 8). all influence the non-linear saturation 



of the MRI in simulations. 

Saturation of the MRI may be related to secondary (par- 
asitic) instabiUties ( Goodman & Xu 1994, hereafter GX94) 
that feed upon the free energy afforded by the MRI 
(see, e.g.. iKnobloch & JulienI 120051: lUmurhan et al.1 120071; 



iTatsuno & Dorlandl 120081: IVishniacf |2009| for alternative 
ideas). These instabilities are often invoked to explain some 
of the behavior observed in numerical simulations. The ab- 
sence of explicit dissipation from most numerical studies and 
from GX94's analysis, however, impedes quantitative inter- 
pretations. 

Both primary and secondary instabilities may be subject 
to non-ideal effects. The growth rates and wavenumbers 
of the fastest growing MRI primaries, and the relative ori- 
entations and magnitudes of their velocity and magnetic 
field pertu rbations, are sensitive to dissipation coefficients 
(see, e.g.,|L esur & Longaretti 20071 ILesaffre & Balbus 2007; 
iPessah & Ch an 2008). Viscosity may slow the growth of the 
Kelvin-Helmholtz secondaries identified by GX94, and non- 
Electronic address: [mpessah@ias.edu, jeremy@astro.princeton.edu | 



zero resistivity may enable resistive instabiliti es such as tear- 
ing modes (see, e.g.. lBoyd & Sandersonl2003h . 

Here we summarize a parametric study of parasitic insta- 
bilities in dissipative regimes accessible to current numeri- 
cal simulations. We adopt the incompressible limit, which is 
relevant for initial fields so weak that saturation occurs with 
sub-equipartition fields. The fastest growing, non-ideal para- 
sitic modes are related to Kelvin-Helmholtz and tearing-mode 
instabilities. They are non-axisymmetric and have horizon- 
tal wavelengths roughly a factor of 2 larger than the vertical 
wavelength of the primary MRI mode. Our findings suggest 
that current simulation domains may bias the saturation of 
the non-ideal MRI b y excluding the fastest MRI and parasitic 
mode s (see also Hawlev et al.lll995; Sano 2007; Bodo et alj 
l2008h . 

2. VISCOUS, RESISTIVE PRIMARY MRI MODES 

Consider an incompressible, Keplerian background with 
constant viscosity v and resistivity 77 and threaded by a ver- 
tical magnetic field, B^, aga inst perturbation s paral lel to the 
background magnetic field. iPessah & ChanI (1200 8') showed 
that in the shearing box, exact unstable solutions exist of the 
form 

:-qno{r-ro)^ + Vosm{Kz)e^\ (1) 
B,z + BoCOs{Kz)e^\ (2) 



V - 
B 



where q = — dlnfi/dlnr and fio is the local Keplerian fre- 
quency. For a given wavenumber K, the growth rate F satis- 
fies the dispersion relation 

{K^vl,+V,T^fW{K^vl,^Tl)-^K^vl,^ll = 0, (3) 



where T,. 



r 



\/2{2 — q)flo is the epicyclic frequency, vaz = Bz/\/4:TTp 
is the Alfven speed, and p is the density. The relative strength 
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Fig. 1 . — Normalized growth rate.s of the fastest growing parasitic modes vs. the orientation of the horizontal wavevector fc^ with respect to the radial (d = 0) 
direction. In each panel, the MRI magnetic field Bq = Bff*-{y, r]) is such that the fastest parasitic growth rate, maximized over fej^, 8 and k^, matches the 
growth rate of the fastest primary MRI mode, rmax(^'i ^), for the indicated ' viscous Elsasser and magnetic Prandtl numbers, A[/ and Prri. The angles 0^ and 
6*3 mark the directions of the horizontal velocity and magnetic fields of the fastest MRI mode. The fastest parasites lie at 8 ^ By or 8 ^ 6^ and are associated 
with Kelvin-Helmholtz and tearing mode instabilities, respectively. Tearing modes gain prominence at lower Pm. 



Vq/Bq of the MRI velocity and magnetic fields, Vq — 
Vo (cos 6'v, sin^v, 0) and Bq = (cos 0b, sin 0b, 0), and 
their directions Oy and ^b, are known functions of {i/, rj, K). 
The growth rate T has a unique maximum, rinax(i^, at 

Numerical simulations impose periodicity lengths on the 
shearing box. We avoid finite-volume effects in our analysis 
by adopting length and time scales based on intensive param- 
eters: Lq = VAz/^o and Tq = l/flo- Viscosity and resistiv- 
ity introduce two new scales, which we subsume into the di- 
mensionless quantities A^, = v\^/vVlo and = v\^/rjflo, 
whose ratio is the magnetic Prandtl number, Pm = A^/A^ = 
i/frj. The quantity Aj, is known as the Elsasser number, while 
its viscous counterpart is related to the Reynolds number 
(see below). Throughout the rest of the Letter, unless other- 
wise mentioned, we use the scales Lq and Tq to define dimen- 
sionless variables. In these units, magnetic field strengths are 
defined relative to the background field B^, while A^ 
and Ari rj^^. 

3. VISCOUS, RESISTIVE PARASITIC MODES 

To facilitate our analysis of the stability of the non-ideal 
MRI modes against parasitic instabilities, we invoke some ap- 
proximations similar to those adopted by GX94 for ideal mag- 
netohydrodynamics (MHD). We assume that the amplitudes 
of the primary modes are large enough that we can neglect 
the influence of the weak vertical background field, the Cori- 
olis force, and the background shear flow on the dynamics of 
the secondary modes. Furthermore, the energy source for sec- 
ondary instabilities should increase with (some power of) the 
amplitudes of the primary modes. Therefore, the secondary 
growth rates should eventually outstrip the primary ones. We 
thus assume, following GX94, that we can neglect the tempo- 
ral variation of the sinusoidal velocity and magnetic fields of 
the primary modes. 

Under these assumptions, the equations governing the dy- 
namics of the secondary instabilities are 



dtSv + (Aw V) 5v + (Sv-V) Av = 
+ {AB-V)6B + {6B-V)AB 
dtSB + (AvV) SB + (SvV) AB = 

(AB-V) Sv + (SB-V) Av + rjV'^SB 



V{5P + AB-SB) 
vV'^Sv , (4) 



where Av = Vq sm{Kz), AB = Bq cos{Kz), V-5v = 
V-(5B = 0, and SP stands for the pressure perturbation. We 
seek solutions of the form 



6v{x, t) = 5vo{z) exp [st — ik-x] , 
6B{x, t) = SBo{z) exp [st — ik-x] 



(6) 
(7) 



where the amplitudes dvo{z) and SBo{z) are periodic in z 
with period 2tt/K. Substituting expressions ^ and (|7]i into 
Equations (|4|i and (|5]), and using the divergenceless nature of 
the perturbed fields, we derive a set of higher order differ- 
ential equations for the vertical components of the secondary 
velocity and magnetic fields. 



{s + vQ)Q5vz-i{k^-Av){Q- 
+ i{k^-AB){Q 
{s + rjQ)SBz + i{k^-AB)Svz 



-K^)SBz = 0, 
- i{kh-Av)SBz 



(8) 
0, (9) 



(5) 



where the horizontal wavevector fch is such that k — k\^ + 
kzZ, and the differential operator Q = — d^. The ratio 
kz/K need not be rational but < k^/K < 1/2 (GX94). 
Equations ^ and with suitable boundary conditions, 
pose an eigenvalue problem' for the parasitic growth rate 
s = s(^, rj, K, Bq, kz, 0, k\i), where 9 denotes the angle be- 
tween flie horizontal wavevector kh and the radial direction. 

iLesur & Longarettil (1200 7*. hereafter LL07) carried out a 
systematic study of the saturation of the MRI in incompress- 
ible MHD with explicit viscosity and resistivity. They pre- 
sented results for the dimensionless stress at saturation for a 
grid of models with Pm = {0.125,0.25, 1.0,4.0,8.0} and 
Re = SLl/iy = {200, 400, 800, 1600, 3200, 6400}, with 
5* = 3r2o/2. To relate our results to those of LL07, we 
translate their Reynolds number Re into A^ = 3Re/2/3. 
Their parameter (3 = S'^L'l/vj^^ is a proxy for the plasma 
f3 parameter in a stratified disk with equivalent height Lz (in 
our units Lz = 2V^/3). Setting /3 = 100, as in LL07, 
it follows that the values for Re cited above correspond to 
Ai, = {3, 6, 12, 24, 48, 96}. Note that these are rather smaU. 

We solved equations ([8]) and (|9]) for all the possible combi- 
nations of the set of values of A^ and Pm defined above and 
searched the parameter space defined by {Bq, kz, 0, fch) in or- 
der to identify the most relevant, fastest growing secondary 

' The method of solution of these differential equations as well as a dis- 
cussion of the physics of the secondary modes will be presented elsewhere. 
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Fig. 2. — Predicted dimensionless stress at saturation, Qsat, as a function of and Pm, if saturation occurs when the fastest parasitic and primary MRI 
growth rates match. The indicated values for correspond to the Reynolds numbers studied by LL07, Rc = {200, 400, 800, 1600, 3200, 6400}. In the left 
panel the fastest MRI and parasitic modes are allowed to evolve unimpeded. In the right panel the primary and secondary modes considered are the fastest MRI 
and parasitic modes that can fit in a domain with (Lr, L^, Lz) = (1, 4, 1) X 2y/j3/3 with /3 = 100, as considered in LL07. 



modes. We have not been able to find unstable parasitic modes 
with fch > K. This result generalizes the findings of GX94 to 
the non-ideal MHD regime. 

4. SATURATION OF THE MAGNETOROTATIONAL INSTABILITY 

The secondary modes will be clearly dynamically impor- 
tant when their growth rates are comparable to, or greater 
than, the growth rates of the primary modes upon which they 
feed. We refer to this instance as the "saturation" of the pri- 
mary MRI mode. It is then convenient to define the saturation 
amplitude ry, A') as the amplitude that the magnetic 

field produced by the MRI must have grown to in order for 
the instantaneous growth rate of the fastest parasitic mode, 
Smax('^, K), to match that of the primary, T{v, rj, K)^. 

In the ideal limit, the growth rate of the secondary modes 
derived from equations (H) and ^ is linear in the amplitude of 
the primary magnetic field (GX94). Thus, the amplitude Bq 
at which the growth rate of the fastest secondary equals the 
growth rate of a given primary mode can be estimated after 
solving these equations with i/ = rj = 0. However, in the 
non-ideal case the amplitude Bq cannot be scaled out of the 
problem, and the growth rate of the secondary modes depends 
on it in a non-trivial way. Figure [T| shows the fastest growth 
rates SmaxC^^, ?7, ^max, ^2 , of various secondary modes that 
feed off the fastest primary MRI mode for three combinations 
of Ajy and Pm, with curves for several values of k^- In all of 
the cases shown, Bq — i?Q^'(z/, 77, i^max), i-C-, the amplitude 
of the primary MRI mode is such that the fastest secondary 
growth rate, s,„ax(i^, Vi K^?>.y.) = '^miix{v, i]). 

For all the cases within the explored dissipative regime, the 
fastest parasitic modes are non-axisymmetric (6 7^ 0), have 
the same vertical periodicity as the primary mode (fc^ = 0), 
and have purely real growth rates. The fastest modes have 
horizontal wavevectors that are nearly aligned with either the 
velocity or the magnetic field of the primary {Av, AB). The 
first type are clearly related to Kelvin-Helmholtz instabili- 
ties while the latter are related to tearing modes. The ra- 
tio between the horizontal wavenumber of the fastest para- 
sitic mode, fch,max(i^, ?/, -K^max), and the wavenumber of the 

^ In neglecting the temporal dependence of the background, it is assumed 
that the secondary growth rates are large compared to the primary growth 
rate; thus our definition of saturation implies an extrapolation to the regime 
where this assumption is not strictly satisfied. 



fastest MRI mode is rather insensitive to either A^, or Pm; 
fch,max/^max varics from 0.59, in ideal MHD, to 0.46, in the 
cases with highest viscosity and resistivity. 

We calculate the dimensionless stress at saturation as 
a,at = f',f/{SL,)\ where f^f = i?-* - M^', is the 
sum of the Reynolds and Maxwell stresses 



psat - 
JX^0 - 



-1 / J/o^,f(z)t/o^,^;(z)dz. 



— / B^^::iz)B^^liz)dz . 



(10) 



(11) 



LJ2 



These expressions are integrated to obtain the dimensionless 
stress ttsat in terms of the parameter /?, 

asat = ^ [(Kf"')' sin 29y - (Bg-')^ sin 20b] ■ (12) 

The left panel of Figure |2] shows asat in the case where 
both primary and secondary instabilities evolve unimpeded. 
There are competitive effects that set the value of asat, which 
is dominated by the Maxwell stress. As dissipation increases, 
the saturation amplitude Bq^' increases from Bg^' — 3.8, in 
ideal MHD, to Bq^^ ~ 5.5, in the cases with high dissipa- 
tion. However, the angle 6*3 decreases toward tt/2, so that 
I sin20B| decreases. Both effects roughly compensate each 
other so that the final value of ctgat changes only by a factor 
of 2. Therefore, for the range of dissipation coefficients that 
we explored, the saturation amplitude of the primary modes 
should be insensitive to the dissipation coefficients if both the 
fastest primary and secondary instabilities are permitted. The 
value of ckgat does not depend on Pm for large A^, and the 
range of Pm for which the results do not depend on Pm seems 
to increase with A^,. 

The right panel of Figure |2] shows ckgat when the pri- 
mary and secondary modes considered are the fastest MRI 
and parasitic modes that can fit in a domain {Lr, L^, L^) ~ 
(1, 4, 1) X 2V^/3 with (3 = 100, as considered in LL07. The 
purpose of this exercise is to emulate the situation in which 
a finite simulation domain might constrain the availability of 
primary and/or secondary modes. The dependence of agat on 
the dimensionless numbers ( A,^ , Pm) is modified with respect 
to the case where there are no limitations on the range of pri- 
mary or secondary modes; agat decreases with decreasing Pm 
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for a wider range of A^. This trend is similar to that observed 
by LL07, but the dependence on Pm is less pronounced and 
the predicted value of Ogat is smaller than that found in non- 
linear simulations (by a factor of 6 at Pm = 1). Most of 
the differences in Ogat associated with domain size are due 
to the limitations on the primaries. Although the fastest sec- 
ondaries do not fit in the domain, there are other secondary 
modes with comparable growth rates that can lead to "satura- 
tion" at slightly larger primary amplitudes. 

5. DISCUSSION 

We have investigated the spectrum of parasitic instabili- 
ties that feed off the MRJ in viscous, resistive MHD, focus- 
ing our attention on the parameter space currently accessible 
to numerical simulations. Our study suggests that important 
differences between two-dimensional and three-dimensional 
simulations are to be expected. The fastest parasitic modes 
are non-axisymmetric and have the same vertical periodic- 
ity as the primary upon which they feed. They tend to have 
wavevectors that are almost aligned with the velocity or mag- 
netic fields generated by the MRJ. As we will detail in a sep- 
arate paper, the fastest parasitic instabilities can be roughly 
grouped into Kelvin-Helmholtz and tearing mode instabili- 
ties. The first type feed off the sinusoidal velocity field of 
the MRI and are quenched when the viscosity is increased. 
The second feed off the MRJ currents and are enabled by re- 
sistivity. 

The values of ctsat that result from the rather crude proce- 
dure that we followed to find the amplitudes of the MRJ fields 
at saturation are similar, within factors of a few, to the values 
obtained from numerical simulations in the turbulent regime. 
Our results suggest, however, that the saturation amplitude of 
the MRI may not be well determined if the simulation domain 
limits the available primary and secondary modes. The do- 
mains should be large enough vertically to accommodate the 
fastest primary mode, i.e., > 2tt / K^a;y^{v, ij), and should 
have aspect ratios that allow the fastest parasitic modes, i.e., 
Lr,L^ > 2L^ (see ©. 

Having gained some insight into the dynamics of primary 
MRI modes and their parasitic instabilities, it is instructive to 
review some recent numerical results: 

• LL07 carried out a series of shearing box simulations 
in incompressible MHD with explicit dissipation. For the 
range of parameters that they explored, 0.1 < Pm < 10, 
the stress at saturation decreases with decreasing magnetic 
Prandtl number A weak dependence on the Reynolds num- 
ber cannot be discounted. For several of the runs, the most 
unstable MRJ mode does not fit within the numerical domain, 
and neither do the fastest parasitic modes. It is worth asking 
whether the observed trends of ctsat with Pm and Re might 
be biased by these constraints. 

• iMasada & Sand (l2008h explored viscous effects in two- 
dimensional simulations. In several cases, "saturation" is not 
achieved since the turbulent stresses are still increasing at the 
end of the runs. Two different effects might be playing a role 
in the observed behavior Axisymmetric parasitic modes grow 
at only a fraction of the rate of the fastest non-axisymmetric 
modes. Kelvin-Helmholtz parasites are further slowed by vis- 
cosity. Therefore, the MRJ field needs to grow to higher am- 
plitudes before the secondary instabilities can compete with 
the most unstable (available) primary mode. 

• iBodo et al.l (12008.) exposed a dependence of saturation on 
the aspect ratio of the simulation domain. They found that the 



saturated stresses decrease when Lr/Lz varies from 1 to 4, 
and found less significant differences between aspect ratios 4 
and 8. We have shown that even in non-ideal MHD, the fastest 
parasitic modes have horizontal wavelengths roughly twice as 
large as the vertical wavelength of the dominant primary MRI 
mode. It is tempting to attribute Bodo et al's (2008) results to 
the exclusion of parasitic modes at the 1 : 1 aspect ratio. 

While this Letter was in draft form, we learned of the in- 
vestigations of Latter et al. (2009, hereafter LLB09), which 
are similar to our own. A brief comparison of methods and 
conclusions is in order Both studies generalize GX94 to re- 
sistive MHD and characterize the types and growth rates of 
the non-ideal parasitic modes. LLB09 neglect viscous effects, 
which we include, but supplement their analysis with direct 
numerical simulations. A major conclusion of both studies is 
that parasitic modes can be limited by the size of the simula- 
tion domain. Thus box size may be as important as numerical 
resolution for the saturation of the MRJ. 

LLB09 ultimately conclude, however, that parasitic modes 
are unimportant for two reasons. First, they suggest that par- 
asites should "overtake" the primaries only when the latter 
reach large non-linear amplitudes, and then the gas behaves 
compressibly. We have shown that the secondary modes be- 
come dynamically important when the MRJ magnetic field 
is of the order of a few times the vertical field. Therefore, 
we argue that if the initial vertical field is sufficiently weak, 
then the parasites can erupt before the magnetic pressure ap- 
proaches the gas pressure, although this regime is expensive to 
simulate with compressible codes. Second, LLB09 find that 
in large boxes — where parasitic modes are possible — channel 
modes emerging from the turbulent regime do not reach large 
amplitudes; they conclude that something other than parasitic 
modes must be responsible for their disruption. This may be 
true, although we suggest another explanation below. While 
LLB09 speculate about multi-mode interactions, they offer 
no definite or calculable alternative to parasitic modes. At 
present, the properties of parasitic modes provide the only an- 
alytical guidance, other than the linear MRJ dispersion rela- 
tion, to the choice of dimensions for incompressible or weak- 
field simulations. 

The numerical simulations seem to suggest that the pri- 
mary modes reach higher amplitudes than predicted in Sec- 
tion |4] Perhaps terms neglected from equations (|8]l and (|9]l 
delay the onset of the secondary instabilities by reducing their 
growth rates. Alternatively, saturation may occur not when 
the secondaries achieve the same growth rate as the primary, 
but when they achieve the same amplitude. The latter occurs 
later when the MRJ develops from a quiet start. Let tg be the 
time at which the growth rates match, when the primary has 
amplitude BQ{tg). The instantaneous parasitic growth rate 
Smax(0 ^ [So(t)/So(ig)]ri„ax for t > tg. Supposc that 
the fastest parasite begins with amplitude e <C 1 relative to 
the primary at tg. Then when secondary and primary am- 
plitudes match, at a later time ta, the primary has grown by 
a factor Bo{ta)/ Bo{tg) ~ 1 + ln(l/e). That is, the equal- 
amplitude criterion predicts a saturation amplitude larger by 
~ ln(e/e) compared to our previous equal-growth-rate cri- 
terion. LLB09 seeded their simulations with small but un- 
specified noise, which set the initial amplitude of the para- 
sites. If e ~ 10"'' then the overshoot factor ln(e/e) « 10. 
In the fully turbulent regime, the parasites would start from a 
larger amplitude, say e ^ 0.1, leading to a smaller overshoot 
ln(e/e) ~ 3. Channel modes are then less prominent in the 
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turbulent state because of this smaller peak amplitude as well 
as tiie larger background of other modes. Further numerical 
investigations should shed Ught on these issues. 

This Letter benefited from trenchant criticism by an anony- 
mous referee. M.E.R is grateful to Chi-kwan Chan, Peter Gol- 



dreich, and Aldo Serenelli for useful discussions and grate- 
fully acknowledges support from the Institute for Advanced 
Study. This work was supported in part by NSF award PHY- 
0821899 "Center for Magnetic Self-Organization in Labora- 
tory and Astrophysical Plasmas". 



REFERENCES 



Balbus, S. A. & Hawley, J. F. 1991, ApJ, 376, 214 
. 1998, Rev. Mod. Phys., 70, 1 

Bodo, G., Mignone, A., Cattaneo, F., Rossi, P., & Ferrari, A. 2008, A&A, 
487, 1 

Boyd, T. J. M. & Sanderson, J. J. 2003, The Physics of Plasmas (Cambridge: 

Cambridge University Press) 
Chandrasekhar, S 1960, Proc. Nat. Acad. Sci., 46, 253 
Fleming, T. P, Stone, J. M., & Hawley, J. F. 2000, ApJ, 530, 464 
Fromang, S. & Papaloizou, J. 2007, A&A, 476, 1113 
Fromang, S. & Papaloizou, J., Lesur, G.,& Heinemann, T. 2007, A&A, 476, 

1123 

Goodman, J. & Xu, G. 1994, ApJ, 432, 213 

Hawley J. P., Gammie, C. F., & Balbus, S. A. 1995, ApJ, 440, 742 

Knobloch, E. & Julien, K. 2005, Phys. Fluids, 17, 094106 

Latter, H. N., Lesaffre, P, & Balbus, S. A. 2009, MNRAS, 394, 715 



Lesaffre, R & Balbus, S. A. 2007, MNRAS, 381, 319 

Lesur, G. & Longaretti, R Y. 2007, MNRAS, 378, 1471 

Masada, Y. & Sano, T. 2008, ApJ, 689, 1234 

Pessah, M. E. & Chan, C. K. 2008, ApJ, 684, 498 

Pessah, M. E., Chan, C. K., & Psaltis, D. 2007, ApJ, 668, L51 

Sano, T. 2007, Ap&SS, 307, 191 

Sano, T. & Inutsuka, S. 1. 2001, ApJ, 561, L179 

Sano, T., Inutsuka, S. 1., Turner, N. J., & Stone, J. M. 2004, ApJ, 605, 321 
Tatsuno T. & Doriand W. 2008, Astron. Nachr., 329, 688 
Velikhov, E. P 1959, JETP, 36, 1398 
Vishniac, E. 2009, ApJ, 696, 1021 

Umurhan, O. M., Menou, K., & Regev, O. 2007, Phys. Rev. Lett., 98, 034501 



